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Far-infrared vibrational spectroscopy by multiple photon dissociation has proven to be a very 
useful technique for the structural fingerprinting of small metal clusters. Contrary to previous 
studies on cationic V, Nb and Ta clusters, measured vibrational spectra of small cationic cobalt 
clusters show a strong dependence on the number of adsorbed Ar probe atoms, which increases 
with decreasing cluster size. Focusing on the series CoJ to Cog' we therefore use density-functional 
theory to analyze the nature of the Ar-Co^ bond and its role for the vibrational spectra. In 
a first step, energetically low- lying isomer structures are identified through first-principles basin- 
hopping sampling runs and their vibrational spectra computed for a varying number of adsorbed 
Ar atoms. A comparison of these fingerprints with the experimental data enables in some cases a 
unique assignment of the cluster structure. Independent of the specific low-lying isomer, we obtain 
a pronounced increase of the Ar binding energy for the smallest cluster sizes, which correlates nicely 
with the observed increased influence of the Ar probe atoms on the IR spectra. Further analysis of 
the electronic structure motivates a simple electrostatic picture that not only explains this binding 
energy trend, but also why the influence of the rare-gas atom is much stronger than in the previously 
studied systems. 



I. INTRODUCTION 

The intricate relationship between structural and elec- 
tronic degrees of freedom of metal clusters in the sub- 
nm size range dictates a precise knowledge of the geo- 
metric structure when aiming to exploit the unique opti- 
cal, magnetic and chemical properties of these materials. 
Experimentally, obtaining information on the atomic ar- 
rangement of such small clusters in the gas phase re- 
mains a challenge. Investigations of the reactivity of 
metal clusters with small molecules, in particular their 
saturation behavior, have been used to draw conclu- 
sions about the clusters atomic arrangements ii Photo- 
electron spectroscopy provides direct information on the 
electronic structure of clusters which can be related to 
their geometries.^ Information on the shape of cluster 
ions can be obtained from ion mobility measurements^!^ 
and, recently, electron diffraction of trapped charged 
clusters has shown to be promising for revealing their 
structure. ^'^'^ Another possibility that is particularly 
sensitive to the internal cluster structure is the measure- 
ment of vibrational frequencies. A corresponding tech- 
nique that has recently been successfully employed to de- 
termine the structure of cationic and neutral metal clus- 
ters containing between three and 20 atoms is far-infrared 
(vibrational) resonance enhanced multiple photon disso- 
ciation (FIR-MPD) spectroscopyi^iiaiiiii^ii^ The essen- 
tial idea of this technique is to irradiate rare-gas com- 
plexes of the targeted clusters with infrared (IR) light 
in the range of the structure-specific vibrational funda- 
mentals. When the IR light is resonant with a vibra- 
tional mode in the cluster complex, the complex can ab- 
sorb several photons and subsequently evaporate off one 
or more rare-gas atoms. Recording the resulting abun- 
dance changes of the rare-gas complexes as a function 
of the IR frequency yields the desired spectra that can 



then be compared to computed IR absorption spectra for 
different isomer structures obtained e.g. using density- 
functional theory (DFT). As long as different isomers 
exhibit distinct vibrational fingerprint patterns this then 
enables a unique determination of the cluster structure. 

A fundamental assumption behind this approach is 
that the rare-gas atoms do not significantly influence 
the vibrational spectrum and merely act as a probe 
for detecting the resonant absorption of IR photons by 
the metal clusterfH A negligible influence of the rare- 
gas atoms was indeed inferred from insignificant differ- 
ences in the IR-spectra of previously studied cationic 
V+ (n = 3 - 23)^, Nb+ (n = 5 - 9>U.iil and Ta+ 
(n = 6 — 20)^ complexes that contained one or more Ar 
atoms. Restricted test calculations comparing the IR- 
spectra of bare and Ar-complexed and clusters 
led to the same conclusionSj^, so that the theoretical 
modeling for the mentioned systems focused exclusively 
on the IR spectra of the bare metal clusters. 

We find the situation to be markedly different for small 
cationic Co clusters. Here, the measured FIR-MPD spec- 
tra show an intriguing dependence on the number of ad- 
sorbed Ar atoms, which is strongest for clusters contain- 
ing less than six cobalt atoms. This prompts us to per- 
form a systematic DFT analysis of the nature of the Ar- 
Co^ bond and its implications for the vibrational spectra. 
For the series Coj" to CoJ we first employ basin-hopping 
(BH)i^ to sample the configurational space and identify 
energetically low-lying isomers. Explicitly accounting for 
the Ar probe atoms, a comparison of the computed IR 
spectra for these isomers with the measured data allows 
in some cases for a unique assignment of the experimen- 
tal cluster structure. The calculations reveal furthermore 
a characteristic increase of the Ar binding energy for the 
smaller clusters, rationalizing the observed increased in- 
fluence of the Ar probe atoms on the IR spectra. We 
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can trace this binding energy trend back to the predom- 
inant contribution to the Ar-Co^ bond arising from the 
polarization of the rare-gas atom in the electrostatic field 
of the cationic cluster. This motivates a simple electro- 
static model that not only reproduces the binding en- 
ergy trend, but also explains why the interaction of Ar 
is much stronger than in the previously studied systems, 
where little influence of the probe atoms on the vibra- 
tional spectra is observed. 



II. EXPERIMENT 

The experiments are carried out in a molecular-beam 
setup that is coupled to a beam line of the Free Elec- 
tron Laser for Infrared experiments (FELIX) at the 
FOM-Institute for Plasma Physics in Nieuwegein, The 
Netherlands.— The molecular-beam machine has been 
presented previousl y ^^i^^'^*^ and only a brief description 
is given here. Metal clusters are produced by pulsed laser 
ablation of a cobalt rod using the 2"^^ harmonic output 
of a Nd:YAG laser and by subsequent condensation of 
the plasma in a mixture of Ar in He. Neutral, anionic, 
and cationic clusters are produced in this process and 
pass through a temperature controllable copper channel. 
Only cationic clusters are investigated in this study. Two 
different experimental conditions have been used. With a 
mixture of 0.1 % Ar in He at liquid-nitrogen-temperature, 
clusters containing six or more cobalt atoms bind one 
and, to a much lesser degree, two argon atoms. Under 
the same conditions small clusters up to Co^ form com- 
plexes with up to five argon ligands. At a higher temper- 
ature of 340 K and using a mixture of 0.3 % Ar in He no 
attachment to bigger cobalt clusters is observed, whereas 
Co^ and Co^ add one and two rare-gas atoms which is 
in line with observations by Minemoto et al}^. 

The molecular beam expands into vacuum and passes 
through a skimmer and an aperture before entering the 
extraction region of a reflectron timc-of-flight mass spec- 
trometer. A counter-propagating pulsed far-IR laser 
beam delivered by FELIX is overlapped with the aperture 
to ensure that all species in the beam that are detected 
in the mass spectrometer have been exposed to the IR 
radiation. The laser pulse consists of a 1 GHz train of ps- 
duration micropulses. The duration of such a macropulse 
is several /iS while its energy ranges typically between 15 
and 30 mJ, depending on the wavelength. When the 
IR radiation is resonant with an IR-allowed transition of 
the cluster complex, sequential absorption of single pho- 
tons can take place. The resulting heating of the cluster 
may induce evaporation of the argon atoms, which leads 
to depletion of the complexes in the beam. IR depletion 
spectra are constructed by recording the ion intensities 
of the metal argon complexes as a function of the FELIX 
frequency; from these, the far-infrared absorption spec- 
tra are obtained as described before. As the detection is 
mass selective, the simultaneous measurement of far-IR 
spectra for different cluster sizes is possible. 



III. THEORY 

The DFT calculations are carried out using the all- 
electron full-potential code FHI-aims^ and the general- 
ized gradient approximation (GGA-PBE)^ to describe 
electronic exchange and correlation (xc). For compari- 
son, some calculations are also performed using the local 
density approximation (PW-LDA)S2,. FHI-aims employs 
basis sets consisting of atom-centered numerical orbitals. 
All calculations reported here are done with the "tier2" 
basis set which contains 67 and 45 basis functions for 
Co and Ar, respectively. As integration grids, we used 
158 and 140 radial shells for Co and Ar, respectively, in 
which the number of integration points is successively de- 
creased from 1202 for the outermost shell to 110 for the 
innermost one. Convergence tests focused on the D2d 
ground-state and the higher-lying isomer of Co^-Ar 
{vide infra). Using the hierarchically constructed larger 
basis sets provided in FHI-aims, as well as denser integra- 
tion grids, we recomputed the relaxed geometric struc- 
ture and vibrational spectra of the two isomers, as well 
as their energetic difference, and found these quantities 
to be converged within 0.01 A, 2 cm~^ and 2 meV, re- 
spectively, which is fully sufficient for the arguments and 
conclusions put forward below. 

Local structural optimization is done using the 
Broyden-Fletcher-Goldfarb-Shanno method^**, relaxing 
all force components to smaller than 10~^eV/A. After 
relaxation, the Hessian matrix is determined numerically 
by finite displacements of all atomic positions by 10~^ A, 
and then diagonalized to get the vibrational modes. The 
corresponding infrared-intensities are obtained by taking 
the derivative of the dipole moments along these modes. 
In order to facilitate the visual comparison to the exper- 
iment, the resulting IR-spectra are finally folded with a 
gaussian line shape function of half- width 1 cm~^. 

In order to identify energetically low-lying isomers we 
tested all structures for the bare clusters that had been 
discussed previously in the literature. For increasing clus- 
ter sizes this approach is likely to miss relevant geome- 
tries. This holds for Co in particular, where one has to ex- 
pect several, differently Jahn- Teller distorted versions of 
the same structural motif. In addition, we therefore ran 
extended basin-hoppingi^ runs to achieve a better sam- 
pling of configurational space. BH runs proceed through 
a sequence of jumps from one potential-energy surface 
(PES) minimum to another. For this, an initially random 
cluster structure is subject to so-called trial moves, which 
correspond to a random structural modification, followed 
by a local geometry optimization. The desired impor- 
tance sampling of energetically lowest-lying isomers and 
the possibility to surmount barriers on multiple-funnel 
type PESs is ensured through a Metropolis-type accep- 
tance rule with which a created trial structure replaces 
the current cluster structure as starting point for the fol- 
lowing trial move. 

As specific implementation we chose so-called single- 
particle trial moves, in which a randomly picked atom is 
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displaced in a random direction by a distance of 3.00 A 
(corresponding to 1.5 times the GGA-PBE computed Co 
dimer bond length). In order to prevent an entropy- 
driven dissociation of the cluster during the BH run, we 
disregarded trial moves as well as local relaxations that 
generate loosely connected or partly dissociated struc- 
tures characterized by an atom having a nearest-neighbor 
distance larger than twice the dimer bond length. Simi- 
larly discarded were moves that place atoms at distances 
of less than 0.5 A from each other. In view of the specific 
electronic configuration of Co, we furthermore explicitly 
included the spin degree of freedom into the sampling 
procedure, i.e. after the trial atomic displacement the 
electronic self-consistency cycle was initiated with a ran- 
dom magnetic moment. For specific isomers we further 
supplemented this by fixed-spin moment calculations en- 
forcing spin-states that had not been found in the BH 
runs. However, in almost all cases these states turned 
out to not be local minima of the spin hypersurface, sug- 
gesting that the chosen spin sampling in the BH runs 
yields the relevant spin minima. 

Equivalent BH runs were performed for the Arm-Co^ 
complexes, randomly displacing both Ar and Co atoms. 
The most favorable structures determined by this com- 
pletely unbiased treatment of the Ar probe atoms were 
always Co clusters with Ar adsorbed at a top site, with 
Ar-Co bond lengths that range from 2.45 A for the C04" 
ground-state to 2.63 A for the most stable isomer of Co^. 



IV. RESULTS AND DISCUSSION 
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FIG. 1: Distribution of cationic cobalt clusters and their Ar 
complexes obtained at source temperatures of 340 K (top) and 
80 K (below). At the higher temperature rare gas complexes 
are only observed up to Co^ , whereas at low temperature also 
larger clusters attach Ar and the small clusters bind multiple 
Ar atoms. The dashed red line in the lower plot connects the 
maxima of the Co^-Ar„ peaks. 



A. Setting the stage 

The recorded mass spectra of the cationic cobalt-argon 
complexes reveal profound differences in the binding of 
the rare-gas ligands to clusters containing more or less 
than six atoms, for both employed complex formation 
conditions (77K/0.1% Ar vs. 340K/0.3% Ar). As shown 
in Fig. [1] at the low-temperature conditions, in which 
the bigger clusters bind only one or at maximum two 
Ar atoms, C04 and Co^ readily attach up to five lig- 
ands. The affinity of these small clusters to the rare-gas 
atoms is so pronounced that they still form complexes 
with one or two Ar atoms even at a temperature as high 
as 340 K. This change of the interaction with the probe 
atoms with cluster size is also reflected in the obtained 
FIR-MPD spectra, which are measured in the range 75- 
350 cm~^ with no IR absorption bands detected at higher 
wavenumbers. The spectra show a clear dependence on 
the number of adsorbed Ar atoms, which is strongest for 
the smallest clusters and is essentially lost for clusters 
containing more than seven Co atoms. 

Compared to linear absorption spectra, FIR-MPD 
spectra are known to be subject to distortions, resulting 
from an incomplete fragmentation of complexes contain- 
ing more rare-gas atoms and (cross) anharmonicities in 
the multiple photon excitation processi^. For the first 



of these effect one has to keep in mind that there is 
no mass selection prior to the fragmentation of a cer- 
tain complex. The spectrum of a given complex thus 
contains both a contribution from the depletion due to 
evaporation of the argon ligand as well as an increase in 
intensity due to incomplete fragmentation of complexes 
containing more argon ligands. This can lead to missing 
bands, especially at the low-frequency end of spectra, 
and is avoided by considering only the spectra of clusters 
with the maximum number of argon ligands. The effect 
of (cross)harmonicities is to change the width and posi- 
tions of existing lines. Both effects have been observed 
in preceding work on cationic V+ (n = 3 — 23)^1^^, Nb+ 
(n = 5 - 9)"'i2 and Ta+ (n = 6 - 20)^^ complexes. How- 
ever, the influence of the rare-gas atoms was typically 
small, with only minor shifts of the frequencies of the ab- 
sorption bands on the order of a few cm~^ and with the 
relative intensities of the observed bands barely affected. 
This situation is markedly different in the case of the 
small cationic Co clusters, where the spectra from dif- 
ferent Ar containing complexes exhibit much more pro- 
nounced frequency shifts and rather strong changes in 
the relative intensities, cf. Figs. [3] and 2] below. This 
finding together with the indicated unusually strong in- 
teraction with the Ar probe atoms prompts us to use 
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FIG. 2: Comparison of computed vibrational spectra of the 
identified most stable isomers of CoJ-Ar to Co^-Ar and using 
the PW-LDA (blue line) or GGA-PBE (red line) functional. 
A common scaling factor is applied, matching the main peak 
of the experimental Co^-Ar spectrum. See text and Figs. |3][7] 
below for a description of the different isomers. 



scaled vibrational spectra of the identified most stable 
isomers of Co^-Ar complexes in the targeted size range 
{vide infra), and compares them to the corresponding 
spectra obtained using the PW-LDA functional. In light 
of the particularly controversial situation for the dimer, 
we did not follow the usual practice to determine the 
scaling factor from the comparison of the computed and 
experimental dimer vibrational frequency, but used the 
main peak of the simple fingerprint spectrum of Co^-Ar, 
cf. Fig. [S] below, instead. However, regardless of which 
reference is actually used, the validity of such a common 
scaling approach would have shown up in basically identi- 
cal spectra for all cluster sizes after scaling the PW-LDA 
and GGA-PBE data, which is obviously not the case. 

In this situation, it makes little sense to strive for a 
quantitative agreement of the absolute frequencies when 
comparing the experimental to the theoretical spectra. 
Additionally considering the mentioned differences be- 
tween linear IR absorption and FIR-MPD spectra, the 
best one can possibly hope for is a match between (rel- 
ative) spacings between different peaks in the spectrum, 
or, more modestly, to identify some characteristic fin- 
gerprint pattern that is unique to a specific isomer and 
which would then allow for a structural assignment. It 
is within this perspective that we augment the detailed 
presentation of the low- lying isomers of CoJ to Co^ iden- 
tified in the BH runs with a comparison of the computed 
and measured vibrational spectra of the corresponding 
Ar complexes. Despite the discussed limitations, we ap- 
ply thereby a global scaling to the theoretical spectra 
simply to enable a better visual comparison. Using the 
main peak of the experimental Co^-Ar spectrum as (to 
some degree arbitrary) reference, this yields a global scal- 
ing factor of 0.881 for GGA-PBE and 0.749 for PW-LDA, 
respectively. 



DFT to specifically investigate the Ar-Co^ interaction 
and its role for the vibrational spectra, focusing in view 
of the experimental data on the size range n = 4 . . . 8. 

In general, there is little reason to expect that DFT 
with present-day local and semi-local xc functionals 
yields an accurate description of the intricate electronic 
structure of such small Co clusters. The difficulty of elec- 
tronic structure calculations to properly account for the 
delicate balance between exchange and correlation, be- 
tween magnetism and chemical bonding, and between 
populating d and s orbitals in these systems is already 
well exemplified by the still controversial predictions of 
the true ground-state of the neutral Co dimer (see e.g. 
ref. [2^ and refs. therein). With respect to vibrational 
frequencies a common procedure to correct for some (sys- 
tematic) errors due to an approximate xc treatment is 
to apply a general scaling factor to the calculated fre- 
quencies to bring the computed spectra in better agree- 
ment with the measured ones. The limitations of this 
approach for the cationic Co clusters is, however, well 
apparent from Fig. [21 which compiles the computed and 



B. Low- lying isomers and vibrational spectra of 
their Ar complexes 

1. Cot 

Our calculations for Co^ reveal as most stable isomer 
structure a distorted tetrahedron with D2d symmetry and 
a magnetic moment of 7 /ie • This still quite symmet- 
ric geometry has two (opposite) bonds of 2.50 A length, 
while the remaining four bonds are 2.17 A. The next low- 
est structure is already 0.32 eV higher in energy, has a 
magnetic moment of 9 fiB and corresponds to an even 
further distorted tetrahedron, in which one bond is elon- 
gated to 2.90 A to yield a kind of butterfly geometry. 
Even further up in energy (0.53 eV) is the lowest-energy 
planar structure, with 9/iB, a ^2h symmetry and bond 
lengths of 2.22 A and 2.38 A. The thus obtained energetic 
order is consistent with the interpretation of photodisso- 
ciation spectroscopy data favoring a three-dimensional 
structure for the cationic Co tetramer— . A correspond- 
ing preference for a pyramidal ground-state was also ob- 
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FIG. 3: Comparison of experimental FIR-MPD spectra (up- 
per panel) with the computed IR-absorption spectra of Ar 
complexes of the identified three energetically lowest isomers 
of Co J: The most stable distorted tetrahedron (second panel 
from top), a butterfly- type tetrahedron (third panel from 
top), and a planar rhombus (bottom panel), see text. Each of 
the theoretical panels additionally displays a schematic repre- 
sentation of the cluster geometry also indicating the position 
of the adsorbed Ar atom in the CoJ-Ar complex, the point 
group symmetry, the total energy difference AE of the bare 
cluster with respect to the most stable isomer, and the Ar 
binding energy Eh in the CoJ-Ar complex. For the lowest- 
energy isomer, we additionally show a representation of the 
cluster geometry indicating the position of both adsorbed Ar 
atoms in the CoJ-Ar2 complex. In the CoJ-Ar4 complex each 
of the four Co atoms is directly coordinated to one Ar atom. 



tained in earlier DFT calculations for neutral C04 clus- 
ters, albeit with a much smaller energetic difference to 
the lowest lying planar isomerj ^^i^^ Theory hence pre- 
dicts similar geometries for the neutral and cationic C04 
ground-state, which is in line with a corresponding inter- 
pretation derived from measured bond energy data^'*. 

From this quite consistent picture of the structural mo- 
tif of the CoJ ground-state it is somewhat unexpected to 
achieve only a rather poor agreement of the computed 
vibrational spectra of its Ar complexes with the corre- 
sponding FIR-MPD measurements. Turning first to the 
experimental data shown in Fig. [31 pronounced differ- 
ences are observed for the spectra derived from C04 - Ar, 
Co4-Ar2 and Co4-Ar5. This concerns frequency shifts 
and changes in the relative intensities of a group of pre- 
sumably three modes in the range 160-210 cm~^, as well 
as the appearance of a strong mode at 110 cm~^ and ten- 



tatively a small mode at 90 cm"-'^ only in the spectrum 
from the complex with five Ar atoms. Such apparent 
missing peaks at the low-frequency end of the spectrum 
from complexes can be rationalized by the a Ar-cluster 
binding energy for the first ligands and the comparably 
low photon energy in the far-IR, leading to an incomplete 
fragmentation^,. The frequency shifts particularly of the 
central peak observed for the group at higher wavenum- 
bers are, however, much stronger than those in preced- 
ing work on cationic clusters from group 5 of the Peri- 
odic Table-ii^iiiii^iii^. In this situation, we note that the 
spectrum of the Co^-Ar5 complex can not be affected 
by an incomplete fragmentation of complexes containing 
even more Ar ligands. From the relative abundances of 
Co^-|-Ar4 and Co^-I-Ars in the mass spectrum shown in 
Fig. [5] we expect that the latter complex is is in fact best 
described as CoJ-Ar4 with an additional rather loosely 
bound Ar atom. Correspondingly, the spectrum derived 
from Co^-Ars should give a most faithful representation 
of the vibrational frequencies of a Co^-Ar4 complex. 

The computed IR absorption spectra of the most sta- 
ble D2d isomer complexed with Ar atoms reproduce the 
strong low-frequency mode at 120 cm^^ quite well, cf. 
Fig. [21 supporting the interpretation that the absence 
of this mode in the measured Co^-|-Ar and Co^-|-Ar2 
spectra results from incomplete fragmentation. Compar- 
ing to the also shown computed spectrum of the bare 
cluster it is furthermore clear that this agreement is only 
achieved because of the explicit consideration of the rare- 
gas atom(s) in the calculations. Only the unexpectedly 
strong binding of the Ar atom of 0.3 eV breaks the sym- 
metry of the bare cluster and leads to the appearance of 
IR-active modes in the low-frequency range at all. This 
unusual strength of the Ar-Co^ interaction prevails for 
up to four ligands, each time filling one of the top sites 
at the four corners of the pyramid. Presumably due to 
the lacking dispersive interactions in the semi-local xc 
functional, we did, however, not achieve to reliably bind 
a fifth Ar atom to the cluster at the GGA-PBE level. 
This supports the interpretation of the measured Co^- 
Ars spectrum in terms of the vibrational frequencies of a 
C04 -Ar4 complex, and we correspondingly show in Fig. [31 
the computed spectrum of the latter for comparison. Un- 
fortunately, the favorable experiment-theory comparison 
for the lowest frequency modes does not carry over to 
the higher- frequency group of bands in the experimental 
data. In contrast to the more or less evenly spaced three 
main peaks in this group, theory predicts a large spacing 
between one lower-energy peak and a higher-energy dou- 
blet. While the computed frequencies correspond quite 
well to the lower and higher frequency main peak of the 
experimental group, either the central and dominant ex- 
perimental peak is completely missing or the splitting of 
the doublet is significantly underestimated in the theo- 
retical spectrum. 

The calculated spectra of the second and third lowest- 
energy isomer also shown in Fig. [31 can neither achieve 
a better agreement with the experimental data, so that 
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a significant population of a higher-energy isomer is not 
readily invoked as reason for the discrepancy. At present 
we are thus unable to provide a clear explanation for 
the discomforting disagreement between measured and 
computed data. Two possible reasons could be either an 
insufficient description provided by the employed GGA- 
PBE functional or strong (cross-)anharmonicities in the 
experimental data, which could both particularly affect 
the splitting of the higher-energy doublet in the D2d spec- 
trum. The latter is very sensitive to the distortions of the 
pyramidal ground-state structure caused by Jahn- Teller 
effects^ and the Ar bonding, both of which might not 
be sufficiently treated at the level of a semi-local xc func- 
tional. Simultaneously, one has to keep in mind, however, 
that the unusually strong Ar bonding in the complex 
that we will further discuss in the next section severely 
enhances the multiple photon excitation aspect of FIR- 
MPD, and makes the comparison to computed linear IR 
absorption spectra at least for the complexes containing 
only one or two Ar atoms less justified. 



2. Cot 

For Co^ we determine as most stable isomers three 
slightly differently distorted versions of a tetragonal pyra- 
mid that are all within a 0.06 eV energy range and possess 
all a magnetic moment of IO/xb- The energetically highest 
of the three has still C2v symmetry, while the other two 
are further Jahn- Teller distorted to Ci symmetry by tilts 
of the apex Co atom away from the symmetry axis of the 
slightly folded rhombohedral base. In all three structures 
there are sets of shorter Co bonds of around ~ 2.2 — 2.3 A 
and more elongated bonds to the apex atom of the order 
of 2.5 A. A more symmetric C4v tetragonal pyramid with 
12^B is found at 0.15 eV above the ground-state, followed 
by the lowest-lying planar structure already at 0.81 eV. 
This strong preference for a distorted three-dimensional 
structure in a high-spin state is again similar to corre- 
sponding findings for the neutral pentamer— i^l, further 
confirming that neutral and cationic Co clusters exhibit 
similar geometriea^^. 

The comparison of the computed IR absorption spec- 
tra to the FIR-MPD data is this time slightly more favor- 
able. Figure [4] shows the corresponding data for one of 
the degenerate distorted tetragonal pyramids (the spec- 
tra of the other two differ only insignificantly), as well 
as for the more symmetric tetragonal pyramid and the 
planar structure. The measured spectra from Co^-Ar 
and Co^-Ars complexes show comparable differences as 
in the case of the tetramer complexes: A peak at the 
low-frequency end of the spectrum is only observed for 
the complex containing more Ar atoms and a group of 
bands in the range 175-230 cm~^ exhibits pronounced 
frequency shifts and intensity changes. From the overall 
shape one would expect at least three components be- 
hind this group and this indeed is what is found in this 
frequency range in the computed spectrum for the lowest- 
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FIG. 4: Comparison of experimental FIR-MPD spectra (up- 
per panel) with the computed IR-absorption spectra of Ar 
complexes of identified energetically lowest isomers of Co^: 
One of the three almost degenerate distorted tetragonal pyra- 
mids (second panel from top), a more symmetric C4v pyramid 
(third panel from top), and a planar structure (bottom panel), 
see text. Each of the theoretical panels additionally displays 
a schematic representation of the cluster geometry also indi- 
cating the position of the adsorbed Ar atom in the CojJ'-Ar 
complex, the point group symmetry, the total energy differ- 
ence AE of the bare cluster with respect to the most stable 
isomer, and the Ar binding energy Eh in the Co^-Ar complex. 
In the Co^-Ars complex each of the five Co atoms is directly 
coordinated to one Ar atom. 



energy distorted pyramids. However, even though only 
containing two main peaks in this frequency range, the 
overall spectrum of the symmetric tetragonal pyramid 
isomer does also fit the experimental data rather well, 
so that we cannot safely discriminate between these two 
(in any case rather similar) structures. What is instead 
quite clear from Fig. [?] is that the planar structure can 
not be reconciled with the measurements. 

Due to the already severely distorted geometry of the 
bare ground-state the influence of the rare-gas atom on 
the computed spectrum is less dramatic than in the case 
of the Co tetramer. Rather than leading to the appear- 
ance of a series of additional IR-active bands, its effect 
is more to shift existing peaks and change their intensi- 
ties. With an again rather strong computed binding en- 
ergy of around 0.3 eV this concerns primarily those modes 
that involve displacements in the direction of the Ar ad- 
sorbate. A particularly pronounced example of such a 
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wavenumber [cm^] 

FIG. 5: Comparison of experimental FIR-MPD spectra (up- 
per panel) with the computed IR-absorption spectra of Ar 
complexes of identified energetically lowest isomers of Co^: 
The tetragonal bipyramid with 15/iB (second panel from top) 
and the capped trigonal bipyramid (third panel from top). 
The spectrum for a second tetragonal bipyramid with 13/ib is 
very similar to the one of the 15/ib pyramid and omitted for 
brevity, see text. Each of the theoretical panels additionally 
displays a schematic representation of the cluster geometry 
also indicating the position of the adsorbed Ar atom in the 
Cog'-Ar complex, the point group symmetry, the total energy 
difference AE of the bare cluster with respect to the most 
stable isomer, and the Ar binding energy Eh in the Cog^-Ar 
complex. For the lowest-energy isomer, we additionally show 
a representation of the cluster geometry indicating the posi- 
tion of both adsorbed Ar atoms in the Coji"-Ar2 complex. 

change in intensity is the computed mode at 260 cm"'^ 
for the spectra from the Co^-Ar and Co^-Ars complexes, 
in nice analogy to the variations observed in the corre- 
sponding experimental spectra. 



The two most stable determined isomers of Cog cor- 
respond to shghtly distorted tetragonal bipyramids, with 
the lowest-energy structure exhibiting a magnetic mo- 
ment of IBfiB and another one 0.10 eV higher with 13/iB- 
Both structures possess Dad symmetry and bond lengths 
around 2.3 A. The third lowest-energy geometry is a 
capped trigonal bipyramid also with 13/iB, but with 
0.56 eV above the ground-state this is already signifi- 
cantly higher in energy. With the DFT calculations for 
Coe by Datta et al^ also arriving at a tetragonal bipyra- 
mid ground-state, the similarity between the structural 
motif of cationic and neutral clusters thus continues also 
for the hexamer. 

The clear energetic preference for the tetragonal 



bipyramid motif is supported by the comparison of com- 
puted and experimental vibrational spectra. The FIR- 
MPD data is shown in Fig. [5] and exhibits this time only 
little differences between the measurements of Co^-Ar 
and C0f}"-Ar2 complexes. In both cases the spectrum is 
simple and consists only of two peaks around 250 cm~^. 
This pattern is fully reproduced in the calculated spectra 
for the two tetragonal bipyramids and certainly inconsis- 
tent with the computed spectrum for the capped trigonal 
bipyramid. In the case of the cationic hexamer we there- 
fore achieve an unambiguous assignment of the structural 
motif. Furthermore, the theoretical spectra for the iden- 
tified ground-state structure either bare or with one or 
two Ar atoms shown in Fig.[5]are virtually identical. This 
is distinctly different to the findings for the smaller Co+ 
clusters, and goes hand in hand with a significant reduc- 
tion of the Ar bond strength, computed as 0.12 eV for 
the 15/iB ground-state tetragonal bipyramid. 



4. Co+ 

In the case of Co^ we identify a capped tetragonal 
bipyramid with Csv symmetry and magnetic moment of 
16/iB as most stable structure. Again, we arrive there- 
with at the same structural motif for the cationic ground- 
state as the preceding work for the neutral clusters^l. 
The bond length in the triangle far side from the apex 
atom is 2.31 A and is thus close to the computed bond 
length of 2.36 A of the uncapped Co^ ground-state struc- 
ture, whereas the triangle near side of the apex atom has 
elongated bonds of 2.43 A. The resulting average bond 
length of 2.34 A for the entire cluster is thus slightly ex- 
tended compared to the corresponding neutral isomer re- 
ported by Datta et al^ with an average bond length 
of 2.29 A. Energetically only insignificantly higher than 
the capped tetragonal bipyramid are two spin-variants 
of a distorted pentagonal bipyramid with C2v symme- 
try, namely one with 14/iB only 0.06 eV higher and one 
with 16/iB only 0.08 eV higher. The structural variety 
is further complemented by the energetically next low- 
est isomer, which corresponds to a capped trigonal prism 
also in a 16/zb high-spin state and 0.15 eV higher than the 
ground-state. Still energetically quite close this is then 
followed by another, slightly differently distorted 16/^b 
Csv tetragonal bipyramid (O.lSeV) and a low-symmetry 
16/iB structure at 0.29 eV. 

With this close energetic spacing it is particularly in- 
teresting to compare to the vibrational spectra from the 
FIR-MPD experiments and see if the there determined 
fingerprint pattern does maybe conform with one of the 
slightly higher-energy isomers. The corresponding mea- 
surements are shown in Fig. [SI exhibiting slightly larger 
differences between the data from the complexes with dif- 
ferent number of Ar atoms than in the case of Co^, yet 
still significantly less than in the case of the smaller clus- 
ters C04 and Co^. The computed IR absorption spec- 
trum for the ground-state tetragonal bipyramid agrees 
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FIG. 6: Comparison of experimental FIR-MPD spectra (up- 
per panel) with the computed IR-absorption spectra of Ar 
complexes of identified energetically lowest isomers of Co^; 
The capped tetragonal bipyramid with IG/ib (second panel 
from top), two spin- variants of a pentagonal bipyramid (third 
and fourth panel from top), as well as a capped trigonal prism 
(bottom panel), see text. Each of the theoretical panels addi- 
tionally displays a schematic representation of the cluster ge- 
ometry also indicating the position of the adsorbed Ar atom 
in the Co^-Ar complex, the point group symmetry, the to- 
tal energy difference AE of the bare cluster with respect to 
the most stable isomer, and the Ar binding energy in the 
Co^'-Ar complex. For the lowest-energy isomer, we addition- 
ally show a representation of the cluster geometry indicating 
the position of both adsorbed Ar atoms in the Co^-Ar2 com- 
plex. 



overall rather well with the experiments, reproducing all 
three higher-frequency bands. The observed influence of 
the Ar atom in the comparison to the bare cluster spec- 
trum is also in line with the experimental trend, i.e. it is 
somewhat larger than in the case of Co^, but less than in 
the case of C04 and Co^. Defering the detailed discus- 
sion to the next section we note that this correlates well 
with the now again slightly larger computed Ar binding 
energy of 0.18 eV. 

A major concern in the comparison to the experimen- 
tal data is the computed strong Ar-induced mode at the 
lower frequency end of the spectrum. For this, one can- 



not resort to the argument of apparent missing peaks 
in the spectra of complexes with one or two Ar atoms, 
as there are no complexes with more Ar atoms stabi- 
lized under the employed beam conditions, cf. Fig. [T] 
Comparing the theoretical spectra of the bare and Ar- 
complexed Csv ground-state isomer it is, however, clear 
that the intensity of this particular mode depends sensi- 
tively on the interaction with the Ar ligand. The second 
most stable adsorption site at this isomer is atop one of 
the basal Co atoms that are on the opposite side from 
the capping atom, and its binding energy is only 60 meV 
lower. In the corresponding spectrum (not shown) the 
mode does hardly show up. Most likely, we would there- 
fore attribute the discrepancy between experimental and 
theoretical data with respect to this mode to subtle short- 
comings of the employed xc functional in describing the 
cluster-ligand interaction. Returning to the question of 
a unique structural assignment, one would probably still 
conclude that the spectrum of the ground-state capped 
tetragonal bipyramid agrees overall best with the experi- 
mental data. Nevertheless, one also has to recognize that 
the fingerprint pattern of the different low-energy isomers 
is not sufficiently distinct to rule out small percentages 
of these isomers in the experiment as well. 



Coi 



The most stable isomer found for Co^ is a double- 
capped distorted tetragonal bipyramid with 17/iB, which 
thus - like already the C07 ground-state - also contains 
the very stable tetragonal bipyramid as a structural sub- 
unit and again coincides with the structural motif de- 
termined as most stable for the corresponding neutral 
cluster^^. The geometry has nearly D2d symmetry, with 
only bond distance differences of the order of 0.01 A 
reducing it to the C2V point group. The next lowest iso- 
mer found is only marginally higher in energy (0.09 eV) 
and corresponds to a double-capped trigonal prism with 
a magnetic moment of 17/iB. This is followed by an- 
other double-capped tetragonal bipyramid of the same 
magnetic moment (0.21 eV), and then several differently 
distorted and energetically virtually degenerate capped 
pentagonal bipyramid structures at around 0.47 eV. 

For this largest studied cluster the comparison of the 
computed and measured vibrational data is unfortu- 
nately again rather inconclusive with respect to a po- 
tential structural assignment. The experimental spectra 
are shown in Fig. [7]and exhibit primarily two main peaks 
with some suggested shoulders at the high-frequency end. 
There are only marginal differences between the mea- 
surements from the Co^-Ar and CoJ-Ar2 complexes, 
which is characteristic for all further FIR-MPD spec- 
tra recorded for larger clusters up to 18 Co atoms (not 
shown). This minor influence of the adsorbed rare-gas 
atom is similarly found in all computed spectra for the 
three most stable Co^ isomers, concomitant with a mod- 
erate Ar binding energy of around 0.15 eV. In all cases 
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FIG. 7: Comparison of experimental FIR-MPD spectra (up- 
per panel) with the computed IR-absorption spectra of Ar 
complexes of identified energetically lowest isomers of Co^: 
The double-capped distorted tetragonal bipyramid with 17pB 
(second panel from top), the double-capped trigonal prism 
(third panel from top), as well as another double-capped 
tetragonal bipyramid (bottom panel), see text. Each of the 
theoretical panels additionally displays a schematic represen- 
tation of the cluster geometry also indicating the position 
of the adsorbed Ar atom in the Cog'-Ar complex, the point 
group symmetry, the total energy difference A_E of the bare 
cluster with respect to the most stable isomer, and the Ar 
binding energy Eh in the Cog'-Ar complex. For the lowest- 
energy isomer, we additionally show a representation of the 
cluster geometry indicating the position of both adsorbed Ar 
atoms in the Co^-Ar2 complex. 



the data obtained is furthermore more or less similarly 
consistent with the measured fingerprint pattern, with 
maybe a slightly better agreement for the ground-state 
isomer. Put more modestly there is at least no com- 
pelling evidence that it is not the latter structure that 
is predominantly observed, and thus there is neither a 
reason to question the energetic order determined at the 
GGA-PBE level or to invoke the kinetic stabilization of 
a significant fraction of higher-energy isomers in the ex- 
periments. Part of the reason for the disappointing in- 
conclusiveness are the rather strong Jahn- Teller defor- 
mations found for the high-spin isomers, which make it 
partly a matter of semantics to distinguish between e.g. 
a double-capped distorted trigonal prism or a double- 
capped distorted tetragonal bipyramid. Without clearly 
distinct symmetries the various geometries exhibit only 



quantitative differences in the IR-spectra, which in turn 
severely limits a technique like FIR-MPD which relies on 
an indirect structure determination through characteris- 
tic fingerprint patterns. 



C. Nature of the Ar-CoJ bonding 

While the comparison of computed and measured vi- 
brational spectra in particular of the smallest studied Co 
clusters leaves a number of open questions and doubts 
with respect to the detailed interpretation of the FIR- 
MPD data, the general trend of the influence of the rare- 
gas atoms with cluster size is reproduced quite well by 
the calculations. In general, the Ar-induced changes of 
the IR-spectra become smaller for larger clusters, cor- 
relating nicely with an overall reduction of the Ar-Co^ 
bond strength from an "unusually" high value of 0.3 eV 
for the lowest lying isomers of CoJ to the more "intu- 
itive" value of around 0.1 eV for the largest studied clus- 
ters. Particularly for the smallest clusters, the complex- 
ation with Ar leads to the appearance of new IR-active 
modes in the spectra, hand in hand with frequency shifts 
and splittings of existing peaks and often an increase in 
the absorption intensity. These effects can be primar- 
ily traced back to Ar-induced symmetry breakings in the 
geometric structure of the cluster, together with an en- 
hanced polarization of the cluster electron density that 
we will further analyze below. Both effects depend obvi- 
ously on the strength of the Ar-Co+ interaction, rational- 
izing the correlation with the computed Ar binding en- 
ergy. The explanation for the observation of a much more 
pronounced dependence of the FIR-MPD spectra on the 
number of Ar atoms compared to the preceding work on 
the cationic group 5 clusters2ii2iiiiii2iiH is therefore cou- 
pled to the understanding of the larger Ar binding energy 
at the smaller Co clusters compared to the e.g. ^ 0.1 eV 
computed also with GGA-PBE for the previously studied 
V1+ cluster; 
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As an important key to understand the particulari- 
ties of the Ar-Co+ binding energy, Fig[5] shows a plot 
of the electron density redistribution when adsorbing an 
Ar atom at the ground-state isomer of C04 . Observable 
is a strong polarization of the Ar atom, which is equiv- 
alently obtained in a corresponding analysis of Ar ad- 
sorption at all larger cluster sizes studied. This suggests 
as the predominant contribution to the Ar-Co^ bond a 
mere polarization of the rare-gas atom in the electro- 
static field of the cationic cluster. Further support for 
this view on the nature of the Ar-Co+ interaction comes 
from a multipole decomposition of the charge density of 
the determined ground-state isomers of the bare Co^ 
clusters. Via Hirshfeld analysis we compute the point 
charge and dipole moment on each atom in the bare Co+ 
cluster, yet at the relaxed geometry of the corresponding 
Co+-Ar complex. In a second step we then compute the 
Hartree-type interaction energy of an Ar atom placed at 
the position in the Co^-Ar complex into the frozen elec- 
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FIG. 8: Electron density difference induced by the adsorption 
of Ar at the ground-state isomer of CoJ, see text. For the op- 
timized geometry of the Co^-Ar complex, the plot shows the 
isosurface corresponding to 0.015 electron accumulation 
(red) and depletion (blue) upon Ar addition, obtained by sub- 
tracting from the full self-consistent electron density of the Ar- 
Co J complex the self-consistent electron densities of the bare 
CoJ cluster and an isolated Ar atom at the corresponding 
atomic positions. For better visualization, the drawn skele- 
ton of black "bonds" at the left of the figure represents the 
cluster geometry, while the black cross marks the position of 
the Ar atom. 



2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 
d(Ar-Co^^) [A] 

FIG. 10: Attractive Hartree-type Ar-Co^ interaction energy 
computed in the electrostatic model, see text. Shown is the 
variation with bond length when moving the Ar atom radi- 
ally outwards in the CoJ-Ar (solid line) and in the Co^-Ar 
(dashed line) ground-state complex. The two dotted vertical 
lines indicate the computed equilibrium bond lengths in the 
Co+-Ar (2.45 A) and in the Co+-Ar (2.63 A) complex. 
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FIG. 9: Ar binding energy at the most stable adsorption 
site offered by the determined ground-state isomers of CoJ 
to Co^ . Compared are the values obtained from the real self- 
consistent calculation (black crosses) and the Hartree-type 
interaction energy obtained when placing Ar into the frozen 
electrostatic field formed by the point charges and dipole mo- 
ments of the bare cluster (red circles), see text. 



trostatic field generated by the computed point charges 
and dipole moments of the bare cluster. The correspond- 
ingly obtained binding energies are compared to the fully 
self-consistently computed Ar binding energies at the dif- 
ferent cluster sizes in Fig. [9l From the rather good agree- 
ment we conclude that the polarization of the Ar atom 
describes indeed the predominant character of the bond, 
with larger multipole moments and the back-reaction of 
the Co charge density on the polarized Ar atom contained 
in the true self-consistent calculation forming only a small 
correction. 



Further disentangling the contributions to the bind- 
ing energy, we actually compute only a rather small at- 
tractive component due to the point charges. With the 
calculated dipole moments on the Co atoms generally 
pointing radially outwards of the cluster, a large attrac- 
tion comes instead from the interaction of the Ar atom 
with the dipole moment on the Co atom to which it is 
directly coordinated in the top site position. This is then 
offset by the mostly repulsive interactions with the less 
favorably oriented dipole moments on all other Co atoms 
in the cluster. These findings hold similarly for all of the 
studied cluster sizes, and even carry over to less stable 
Ar adsorption sites and other low-energy isomers than 
the ground-state, for all of which we obtain a similarly 
good agreement of the Ar binding energy computed in the 
electrostatic model with the true self-consistent value. 

On the basis of the thus established electrostatic pic- 
ture we achieve a rather consistent explanation for the 
observed trend of the Ar binding energy with Co cluster 
size, as well as for the stronger interaction compared to 
the previously studied cationic clusters from group 5 of 
the Periodic Table. Taking the ground-state isomer of 
CoJ and CoJ as examples for a stronger and weaker Ar 
bonding, respectively, Fig. [TOl shows the variation of the 
interaction energy in the electrostatic model when chang- 
ing the Ar-Co^ bond length by moving the Ar atom ra- 
dially away from the cluster. Immediately apparent is 
the smaller attraction in the case of Co^, which one can 
loosely rationalize by the larger number of unfavorably 
oriented dipole moments at the other Co atoms of the 
larger cluster to which the Ar atom is not directly coor- 
dinated to. Neglecting electronic self-consistency effects, 
the equilibrium Ar-Co^ bond length is determined by 
the interplay of this attractive electrostatic contribution 
and a strongly repulsive short-range component arising 
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from Pauli repulsion. Since the latter is primarily due 
to the inner shell electrons of the directly coordinated 
Co atom, we can assume this component to be rather 
similar for Co^^ and Co^ to a first approximation. With 
the stronger electrostatic attraction this then leads to a 
shorter Ar-Co^ bond length than in Ar-Co^, and in turn 
to a stronger Ar binding energy at the smaller cluster. 

In the same spirit one can also attribute the compara- 
bly stronger interaction of Ar with Co^ than with the 
previously studied V+, Nb+ and Ta+ clusters to this 
interplay of electrostatic attraction and Pauli repulsion. 
The distance at which the latter sets in depends of course 
on the size of the atomic radius, which is smaller for Co 
than for any of the group 5 elements. With the steep 
Pauli repulsion setting in at already larger distances, one 
can in general expect longer equilibrium Ar-metal clus- 
ter bond lengths and therewith smaller Ar binding ener- 
gies for the group 5 elements. As suggested by Fig. [TUl 
the effect of unfavorably oriented dipole moments is also 
smaller at such increased interaction distances, so that 
one would furthermore conclude on a smaller variation 
of the Ar binding energy with cluster size as in the case 
of Co. Simplistically equating the Ar bond strength with 
the influence of the rare-gas atom on the measured FIR- 
MPD spectra as discussed above, the electrostatic picture 
can thus fully account for the general trends observed in 
the corresponding experiments on these materials. Tac- 
itly assuming a similar character of the rare-gas— metal 
cluster bond, one would even predict that equally pro- 
nounced effects of the probe atoms as in the case of 
Co will be obtained whenever studying cationic clusters 
formed of other elements with small atomic radii or for 
more polarizable heavier rare-gas atoms. 

V. SUMMARY AND CONCLUSIONS 

Summarizing we have used density-functional theory 
to study energetically low-lying isomers of cationic Co 
clusters containing from four to eight atoms. Sup- 
plementing our search for structural motifs with flrst- 
principles basin hopping sampling the identified ground- 
state structures agree for each cluster size with the geom- 
etry determined in preceding DFT studies for the corre- 
sponding neutrals. All structures are in high-spin states 
and exhibit sizable Jahn- Teller distortions. Our specific 
motivation has been a detailed comparison of the vibra- 
tional properties of the clusters with experimental far- 
infrared multiple photon dissociation data. While this 
comparison allows only in few cases for a unique struc- 
tural assignment and leaves some open questions and 
doubts in particular for C04 , there is at least no com- 
pelling evidence that would question the energetic order- 
ing of the isomers obtained at the GGA-PBE level or re- 
ciprocally indicate a kinetic stabilization of a significant 
fraction of higher-energy isomers in the experiments. 



Reproduced quite well by the calculations is the gen- 
eral trend of the influence of the rare-gas probe atoms 
used in the measurements, with Ar-induced changes of 
the IR-spectra becoming smaller for larger clusters. This 
correlates nicely with an overall reduction of the calcu- 
lated Ar-Co+ bond strength from 0.3 eV for the lowest 
lying isomers of Co J to around 0.1 eV for the largest 
studied clusters. This trend can be well explained by the 
nature of the Ar-Co^ interaction, which is predominantly 
characterized by a polarization of the rare-gas atom in 
the electrostatic fleld of the charged cluster. Within the 
corresponding picture we can also rationalize the less pro- 
nounced influence of the probe atoms in preceding FIR- 
MPD measurements on small cationic clusters formed of 
group 5 elements, identifying the smaller atomic radius 
of Co as an important factor. 

For the investigated Co^ vibrational spectra, com- 
plexation with Ar leads in general to the appearance of 
new IR-active modes, hand in hand with frequency shifts 
and splittings of existing peaks and often an increase in 
the absorption intensity. These effects can be primarily 
traced back to Ar-induced symmetry breakings in the ge- 
ometric structure of the cluster. The latter become obvi- 
ously stronger with increasing Ar binding energy, but can 
never be excluded. For a cluster geometry on the verge of 
a Jahn- Teller distortion even the tiniest disturbance by 
a rare-gas ligand may induce symmetry breakings that 
can then substantially affect the spectrum. In this re- 
spect, the present work casts some doubts on the hith- 
erto employed practice of comparing experimental data 
to computed IR-spectra of bare clusters. The influence 
of the probe atom can not be judged from a comparison 
of experimental spectra from different complexes, as the 
latter reveal only the differences between adsorption of 
one or more rare-gas atoms, and not with respect to the 
bare cluster. A computed negligible effect of Ar on the 
spectrum of one cluster isomer does neither justify to dis- 
miss the ligand in the calculations, since this tells nothing 
about the liability towards symmetry breaking of other 
isomers of the same or other cluster sizes. In this situa- 
tion, explicit consideration of the rare-gas atoms in the 
modelling is always advisable at little additional compu- 
tational cost. 
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